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Abstract 

I propose the Langevin equation for 3-geometries in the Ashtekar's formahsm to de- 
scribe 4D Euchdean quantum gravity, in the sense that the corresponding Fokker-Planck 
hamiltonian recovers the hamiltonian in 4D quantum gravity exactly. The stochastic time 
corresponds to the Euchdean time in the gauge, = 1 and A^* = 0. In this approach, 
the time evolution in 4D quantum gravity is understood as a stochastic process where 
the quantum fluctuation of "triad " is characterized by the curvature at the one unit time 
step before. The lattice regularization of 4D quantum gravity is presented in this context. 



^On leave of absence from Department of Physics, Faculty of Science, Shimane University, Matsue 
690, Japan, e-mail address: nakazawa@nbivax.nbi.dk, nakazawa@psl.yukawa.kyoto-u.ac.jp 



In the course to study QCD^ in terms of the Nicoh-Langevin maps it is shown 

that stochastic quantization of 3D Chern-Simons theory recovers the time evolution 
in 4D Euchdean Yang-Mills theory where the stochastic ( or fictitious ) time is inter- 
preted as the Euclidean time 0. There is also another example, found in stochastic 
quantization of matrix models with loop space hamiltonian at the double scaling limit, 
that the stochastic time corresponds to the Euclidean time coordinate in 2D quantum 
gravity 00. These facts motivate us to interpret the stochastic time as the Euclidean 
time coordinate even in 4D quantum gravity. In this short note, I point out that the time 
evolution in 4D quantum gravity is described by a Langevin equation for 3-geometries in 
terms of the Ashtekar's canonical field variables by showing that the corresponding 
Fokker-Planck hamiltonian operator recovers the hamiltonian of 4D Euclidean quantum 
gravity exactly |p. The stochastic time is interpreted as the Euclidean time with the 
gauge = 1 and A^* = 0, where and A^* are the lapse function and the shift vec- 
tor, respectively. The Hartle-Hawking type boundary condition is naturally imposed in 
this scheme by specifying the initial probability distribution functional. Then I propose 
the lattice regularization of this approach with 3D cubic lattice. The corresponding 
F-P hamiltonian defines the lattice regularization for the hamiltonian of 4D Euclidean 
quantum gravity in Ashtekar's formalism. 

At first, I introduce the basic Langevin equation for 3-geometries in the Ashtekar's 
variables to recover the hamiltonian of 4D quantum gravity with the corresponding F-P 
hamiltonian defined latter. The simplest form of the Langevin equation is defined by, 

A4«(x,t) = ACr(x,r) , 
<AC(a:,r)ACj(2/,r)>c = ^/^re'^'^ < F^^^x^r) >^ 5\x ~ y) , (1) 

where A1{x^t) is a SU(2) gauge field in Euclidean Ashtekar's canonical formalism^. In 
this note, Latin indices "i,j,k,..." denote the spatial part of the spacetime coordinate in- 
dices. While the Latin letters "a,b,c,..." denote the spatial part of the internal indices. 

1 I use the notation, = ^^A°; - d.Af + e"'"' A'' A", , and Df" = S^'di + e'^^'^A^, . The field variable 
Af{x) is real in the Euclidean Ashtekar's formalism. 



1 



x,y, ... denote spatial spacetime coordinates. The one step time evolution is defined by 
AA^{x, t) = r + At) — r) in (1). At is the unit of the discretized stochastic 
time and r denotes the stochastic time after n steps, r = uAt. The discretization of the 
stochastic time is considered for convenience of stochastic calculus and for understand- 
ing the precise meaning of the noise correlation. The coupling constant k is defined by 
K = IGnG with G, the gravitational ( Newton's ) constant in the natural unit h = c = 1. 
The noise variable in (1) is not a simple white noise. The expectation value of the R.H.S. 
of the noise correlation is understood to be taken with respect to the noises up to the 



one unit time step before, r — At, in the sense of Ito's stochastic calculus [10]. More 
precisely, if the expectation value is taken only at a specified stochastic time, r, then the 
R.H.S. of the noise correlation is not the expectation value since it does not depend on 
the noise at r. It is also equivalent to require < AQ{t) >= in Ito's calculus. 

The basic Langevin equation is manifestly covariant under the SU(2) local gauge 
transformation, Af{x) ^ A'^ + Df^uj^{x) . While it is not covariant under the spatial 
general coordinate transformation, A°;{x) — > ^4^ + F^^^^{x) , due to the divergent term, 
KAr5^(0)Fj^j^-' , which appears in the transformation of AA^. This divergent term is 
formally ( but not well-defined ) cancelled by adding a term in the R.H.S. of the Langevin 
equation (1). The extra term comes from the invariant path-integral measure as we will 
see later. 

In terms of the solution of the Langevin equation, the following equality holds. 

<Ii,,,a6{At ^{x,t) - Af''''''\x))>^ 

= ^ J^{iriitial)\^^-THpp[A,-k]\^j^{final) 

In the L.H.S., j4^^(a;,r) denotes the solution of the Langevin equation with the initial 
condition, A^(x,0) = Af'^''"°'^\x) . In the R.H.S., Hfp[t^,A\ is defined by. 



HFp[k,A] = -'^Jd'xe^'^K{x)7il{x)F^^{x) 



K 



HMA^] = -i jd'xe^'^mx)K{x)7^l{x) . (3) 



To show the equahty (2), the commutation relation, 

,A){y)] = 5',S'^5\x-y) , (4) 

and the vacuum, |0 > with TT^IO >=< 01^5 = 0, are assumed. The representation of the 
states is given by, 

<A\ = < 0| exp[ J d^xA'l{x)7ri{x)] , 

\A> = U^,aA^n^) - AKx))\0 > . (5) 

The Fokker-Planck hamiltonians (3) are just the hamiltonian for 4D quantum gravity 
without cosmological term in Ashtekar's variables with different operator orderingsQ. 
The annihilation operator tt* is interpreted as the "triad "in the Euclidean Ashtekar's 
formalism. NamelyQ, 7r^(x) = —^e{{x) . The stochastic time evolution with the Langevin 
equation (1) corresponds to the Euclidean time evolution in 4D quantum gravity with 
the gauge fixing, = 1 and A^* = 00. In the stochastic process (1), the "triad "plays 
the role of the noise variable and the "quantum fluctuation "of the triad is character- 
ized by the curvature. Though the two expressions in (2) are precisely equivalent, the 
second definition of the F-P hamiltonian in (3) has an advantage for further discus- 
sion. This is because the vacuum satisfies the ( local ) hamiltonian constraint with 
the operator ordering in Hpp, < 0\H{A{x),7r{x)) = Ti.{A{x),7r{x))\0 >= , where 
Hfp[A,7i] = Jd^xT-C{A{x),7r{x)) . It should be noticed that the operator ordering of 
the F-P hamiltonian operator Hpp is deferent from that appeared in the F-P equation 
for the probability distribution. 

Let us consider the initial distribution dependence of the probability distribution 
functional by averaging the expectation value (2) with respect to the initial probability 
distribution. It is defined by integrating out the initial configuration A,"^*™**"'-* (x) , on 
^ For a detailed derivation of the F-P hamiltonian operator from the Langevin equation, see also 
Ref. §. 

ej, = (?^^^e^ , where q — det{qij) with spatial metric g'-' — e^e^ . 
^ In 2D quantum gravity, it is observed that this gauge fixing recovers the time evolution defined in 

non-critical string field theories |Tl|]. 
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which the solution of the Langevin equation v4"^(x, r) depends, with the distribution, 
p[^(imtiaO 0], in the L.H.S. of (2). The integration of Af'"-'^'"'-^ gives a generalized form of 
the distribution functional P[A, r], which is defined by < 0[A^](r) >^= JVAO[A]P[A, r], 
as follows 

p[A, t] = J PA''^^*^"' < A'^^'^'^^IPIA, 0]e-^^^^[^'*l |A > . (6) 

For an arbitrary observable 0[A], the average with respect to the initial value distribution 
also gives, 

< 0[A^{r)] yp^initi"' < A»"^*^»'|p[i,0]e-^-^^^[^-*lO[i]|0 > . (7) 

In the definition of the expectation value in the L.H.S. of (7), the average is also taken 
with respect to the initial values with the distribution p[A*'*«*««'^ o]. For example, eq.(2) 
is given for 0[A] = U^^a,iS[Af{x) - ^■^•'^''"''^(a;)), with the initial distribution, P[A,0] = 

Il^..,j(Anx)-Af''''''''Hx 



Eq.(7) leads the time evolution equation for the expectation values of observables, 
d S 

- < 0[A^{x,t)] >c= - < H^A^x), ^^]0[A"(^)]U^(.)=A.^(.,.) >c , (8) 

as it is expected from the Langevin equation (1). I here notice that the initial condition 
dependence in (7) implies a constraint for initial distribution itself due to the equation, 

^ < 0[A^{x,t)] -ypA*"P[A^",0]i^i.p[A*",^] < A^'^|e-^^^^[^'*l0[i]|A > . (9) 

The existence of the equilibrium limit requires that the R.H.S. in (9) should be zero at the 
infinite stochastic time. A trivial solution of this constraint is the vanishing curvature at 
any spatial points, -P[A, 0] = Ilx,a,ijS(^Ftj{^)) ■ This initial distribution, however, should 
be excluded because the R.H.S. of (8) is identically zero even at finite stochastic time 
and there is no time development. In general, if we choose a solution of the hamiltonian 
constraint as the initial value distribution, obviously there is no time evolution in the 
Fokker-Planck type equation (8). Especially, the constraint does not allow us to solve 
the Langevin equation with the initial condition A^(*"'**"^) (^x) — for the existence of the 
time evolution. 



4 



To specify the physical boundary condition and to study a class of solutions for the 
hamiltonian constraint of 4D quantum gravity in this context, we may choose the initial 
condition which generates a nontrivial time evolution 

Ph-h[A, 0] = U,^,,,a,iS{Ft^{x))) . (10) 

The spacetime point, x = at r = 0, is identified to the point where the 3D spatial 
manifold is absorbed into nothing in an analogous sense of the Hartle- Hawking type 
boundary condition |T^. 

As it is clear from (7), initial distributions and observables are the key quantities to 
specify 3-boundaries on quantized 4D spacetime manifolds in this approach. These quan- 
tities are the solution of the momentum constraint and the SU(2) Gauss law constraint. 
There are some candidates which may be useful to characterize 3-boundaries, such as 
the extrinsic curvature term, 3D Chern-Simons term, topological invariants [1^ and loop 
variables for which the present formalism can be applied. Under the appropriate 



choice of these observables, one can show the gauge invariance ( gauge independence ) 
of the expectation value of observables, following the standard method to fix the gauge 



in the Langevin equation |[T^[0. For spatial general coordinate invariance and SU(2) 



gauge invariance, the gauge fixed Langevin equation is given by, 

AA^ix, r) = -AriDf^" + Ft^<^^}{x, r) + AQ{x, r) , (11) 

with the same noise correlation as eq.(l). The gauge fixing corresponds to choose the 
multiplier fields, '^"•{x) and $*(x), as specified functionals of A"'s. Then, the F-P hamil- 
tonians are given by, 

Hfj!-\n, A] = Hfp + J £x[<l>\x)Ct\x) + ^'^(x)Q(x)} , 

Hi,%^-\A,7r] = Hfp- Jd^x{<l>\x)C!'\x) + ^%x)C^{x)} , (12) 



instead of eq.(3), with 



Q(x) ^ Dfirl. (13) 



C*^ and Cq defines the momentum constraint and the Gauss law constraint, respectively. 
In the first expression in (12), C**^ and Cq, are defined by replacing the operator orderings 
in the original constraints so that the conjugate momentum, tt^'s, are gathered to the 
left to v4"'s in the constraints. The probability distribution and the expectation value 
defined by the Langevin equation (11) are given by eqs.(6) (7) together with the F-P 
hamiltonian (12) instead of (3). As it is clear from (7) with the second expression of the 
F-P hamiltonian in (12), the insertion of constraints does not change the expectation value 
of the observables provided that the observables are the solution of these constraints, such 
as the extrinsic curvature term and 3D Chern-Simons term. We notice that the operator 
ordering in the second expression in (12) is essential for the local gauge invariance in our 
definition of the vacuum, namely the vacuum is gauge invariant, < Q\C^ ^ = C^ ^lO >=< 
Ol^G ~ ^g\^ ' with the operator ordering in (13). Therefore, (11) gives the 
same expectation value for arbitrary observables as one given by (1), provided that these 
observables are the solution of the constraints in (13). The gauge fixing in (1) corresponds 
to = 1, = = 0, where is essentially the shift vector. On the other hand, as for 
the partition function defined by the Langevin equation (11), which is given by (6) with 
(12), the gauge fixing procedure leads a gauge invariant ( more precisely BRS invariant ) 
partition function with a suitable choice of multiplier field, and as functionals of 
dynamical variables. To fix the gauge more generally for the reparametrization invariance 
of the time coordinate, one has to introduce the lapse function by multiplying the R.H.S. 
of the noise correlation (1) with it and choose it as a specified functional of A^'s. Though 
the problem is not addressed in this note explicitly, it is possible following the general 
framework ||16 |. 

Now let us consider the regularization procedure in this context. As it has been 
noticed in (1), an extra term is necessary for the covariance of the Langevin equation 
under the spatial general coordinate transformation. It is a direct consequence of the 
fact that the Ito's calculus picks up the Jacobian factor which comes from the change 
of variables in the path- integral measure [|1^[^. One way to specify the path- integral 
measure in configuration space is to introduce a regularization for the noise correlation 
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in (1). I here propose a lattice regularization of the Langevin equation and the noise 
correlation (1) in which the invariant property of the regularized Langevin equation in 
the sense of Ito's calculus naturally introduces the invariant path-integral measure in the 
configuration space of link variables. 

The lattice regularized Langevin equation is defined by 

aui{t) = {AWi{T) + AT\g\'/'{DLf,(\g\-'/'g;^^)T^}Ui{^ 
\g\ = det{g^^). 

(14) 

Here the dynamical variable and the noise variable, Ui and AWi respectively, have 
been assigned on the link, / = {x,i), of 3-dimensional cubic lattice, which is specified 
by the site x and its nearest neighbor in the i-th direction x + i. Ui is an element of 
SU(2) group in the adjoint representation. While the noise, AWi is algebra-valued in 
the adjoint representation, AWi = AW^fT", with [T" ,T^] = e^^'^T^ {Dl)i is the left Lie 
derivative on the group manifold defined by 

[{Dr.)t, U,] = S^'T'^Ui , [{Dr.)t, {D^fi] = -e^'\Dr)1 . (15) 

The regularized Langevin equation describes the one step time evolution of the link 
variables, Ui{t + At) = Ui{r) + AUi{t) . It should be noticed that Ui(t + At) is also an 
element of SU(2) group. The quantity gfll , is interpreted as the inverse of the "superspace 
"metric, here the superspace is spanned by the link variables {Ui}. The inverse of the 
superspace metric is given in the following regularized correlation of the noises defined 
on the links, I — (x, i) and /' = {U)])- 

< {AWi)^p{t){AWi,),s{t) >= At < g^^ {t)T:pT'^s > , (16) 

where, 

i^i^jU)^^ = ^{{U{x,i)U{x + l,j)U{x + j,i)^U{x,j)^)^^ 

- {U{x,j)U{x + j,i)U{x + ljyU{x,i)^)^^} . (17) 
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go is an arbitrary constant. 

The lattice regularization has been determined from the following requirements, (i) 
The Langevin equation (14) and the noise correlation (16) should be covariant under the 
SU(2) local gauge transformation, 

U{x,i) ^ V-\x)U{x,i)V{x + i) , 
AW{x,i) ^ V-\x)AW{x,i)V{x) . (18) 

(ii) The model recovers the Langevin equation and the noise correlation in (1) at the naive 
continuum limit defied by, U{x, i) ^ 1 + A\T^a and AW{x, i) ^ AC-T''a. At the limit, 
the coupling constant and the stochastic time are scaled, goo^ I and At —>■ a^^Ar, as 
the lattice spacing a goes to zero, a — > 0. (iii) The regularized Langevin equation (14) is 
covariant under the "general coordinate transformation "in superspace, {U{x,i)}. 

U{x,l)^p^V{x,t)af3[U] , (19) 

where V{x,i) is also an element of SU(2) group in the adjoint representation and an 
arbitrary functional of U{y,j). The second term in the R.H.S. in (14) is necessary for 



the invariance in Ito's calculus |T^. The requirement (iii) is introduced to realize the 
spatial coordinate invariance in the lattice regularization. By using these requirements, 
it is straightforward to construct the model in the fundamental representation. It is given 
by simply replacing the inverse of the superspace metric (17) to the following one. 

Owi'r)TaA = 9o^x,y{Sl3y{ni^jU)a5 " <^a5 (Hi^j [7)^/3 } . (20) 



By applying the argument in Ref. |T^, I have introduced the second term in the R.H.S. 
of (14) for the covariance of the lattice regularized Langevin equation under the trans- 
formation (19). The similar argument is formally possible even in the non-regularized 
version (1) if we identify the R.H.S. of the noise correlation in (1) to the inverse of the 
continuum superspace metric. Indeed, from the requirement (ii), the naive continuum 
limit of these equations, (14) (15) and (16), coincides with those in eq.(l) except a term 
which comes from the second term in the R.H.S. of the Langevin equation (14). The ex- 
tra term has the same structure as the second term in the regularized Langevin equation 
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(14) if we identify the noise correlation in (1) as the inverse of the superspace metric in 
Ashtekar's configuration variables. The extra contribution which comes from this term 
under the spatial general coordinate transformation formally cancels the divergent term, 
kAt6^{0)F^^C,-' , which appears in the transformation in AA'}. The extra term appeared 
in the naive continuum limit of (14) actually represents the contribution from the path- 
integral measure and it is necessary for the invariant property of the Langevin equation 
(1) in a formal sense, however, it is not well-defined without regularization. A more de- 
tailed discussion, especially on the covariant property of the lattice regularization under 
the spatial general coordinate transformation in the present approach, will be published 



in elsewhere UTS 



By using the same method developed in the continuum version, one can derive the 
regularized version of the expectation value, such as (6). I only comment on the three 
important consequences of the lattice regularized Langevin equation (14). The first is 
the corresponding lattice regularized F-P hamiltonian which is given by, 

Hreg = -^Ei^r^'(^^)r{i^r^/'^.i'™f.} . (21) 
^ 1,1' 

The F-P hamiltonian defines the lattice regularization of the hamiltonian for 4D Eu- 
clidean quantum gravity in Ashtekar's formalism^. The second is the invariant measure 
which appears in the path-integral representation of the regularized expectation value of 
an observable (6) after integrating out the canonical momentum variables, (-Dl)^. 

VU = g-^'^IiidUi , (22) 

where dU denotes the left invariant Haar measure. The measure is invariant under the 
general coordinate transformation in superspace (19). 

The third is the Schwinger-Dyson equation in this context. It is given by 

< \g\^l^[{DL)\,{\Q\-^'^Qtv)}T''Ui >= , (23) 



^ In Refs. [|l^|20|, similar lattice regularized hamiltonians have been discussed without the contribu- 
tion of \Q\. 
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at an equilibrium limit, limT-_^oo < U >= 0. 

In this note, I have pointed out that there exist a Langevin system for 3-geometries 
which describes the time evolution in 4D Euclidean quantum gravity. There have been 
some attempts to define 4D quantum gravity in the context of stochastic quantiza- 
tion lH^f^f^l^ll by introducing the stochastic time as the fifth time coordinate, how- 
ever, the philosophy of the present approach is different from these previous works. The 
strategy I would like to adopt here is to characterize the 3-boundaries in the 4D space- 
time manifold by using the solution of both the momentum constraint and the Gauss 
law constraint to prepare the initial distribution and the observables. Then the remain- 
ing constraint is only the hamiltonian constraint which may be realized at the infinite 
stochastic time limit provided there exists an equilibrium limit. This is actually the idea 
which we have learned in stochastic quantization of matrix models in loop space for 2D 
quantum gravity, where the stochastic time is the proper Euclidean time and the hamil- 
tonian constraint is realized as the Virasoro constraint in loop space In order to 
work in this program, one has to construct observables in the lattice regularization. It 
should be noticed that the Wilson loop is not a spatially general coordinate invariant 
observable in the present lattice regularization WE . 



In the Langevin equation (1) and (14), the "triad "in Ashtekar's variables are realized 
as noise variables, which presumably represents the stochastic 3-geometries. The lattice 
regularized Langevin equation (14) may provide a possible basis for numerical simulation 
in 4D quantum gravity. The problem of the present scheme is that the noise correlation 
in the basic Langevin equations, (1) and (14), are not positive definite. It would force 
the Langevin equations to be complex. The field variables, though they are real in the 
Euclidean Ashtekar's formalism, would also become complex in the time evolution. The 
point would be main difficulty for the numerical analysis in this scheme. One way to 
deal with the problem may be to extend the Langevin equation (1) to a class of more 
general gauge fixing. It is always possible by multiplying the noise correlation in (1) with 
lapse function. Then one may chose the lapse function so that the noise correlation keeps 
the value to be positive definite. It is an open question if this gauge fixing procedure. 
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a choice of non-trivial lapse function, make sense in the lattice regularization. It is also 
an open question how we introduce the cosmo logical term in (1) and (14). Apart from 
these questions, the description with the Langevin equation has a topological feature 
in the sense of Nicoli-Langevin map. Such a topological feature would relax another 
difficulty, renormalizability of quantum gravity. I hope that the approach is useful for 
deeper understanding of quantum gravity. 



Acknowledgements 

The author would like to thank J. Ambjorn, J.Greensite, A. Krasnitz, H. B. Nielsen 
and M. Weis for enlightening discussion and comments. He also wishes to thank J. L. 
Petersen for continuous encouragement and all members in high energy group at Niels 
Bohr Institute for hospitality. 



References 

[1] G. Parisi and N. Sourlas, Phys. Rev. Lett. 43(1979)744. 

[2] H. Nicolai, Phys. Lett. B89(1980)341; Nucl. Phys. B176(1980)419. 

[3] G. Parisi and Y. Wu, Sci. Sin. 24(1981)483. 

[4] M. Claudson and M.B.Halpern, Phys. Lett. B151(1985)281. 

[5] A. Jevicki and Rodrigues, Nucl. Phys. B421(1994)278 



[6] N. Nakazawa, NBI-HEP-94-54, |hep-th/941123^ . Modern Phys. Lett. A, in press 



[7] A. Ashtekar, Phys. Rev. Lett. 57(1986)2244; Phys. Rev. 36(1987)1587. 

T. Jacobson and L. Smolin, Phys. Lett. B196(1987)39. 

[9] N. Nakazawa, NBI-HEP-95-26, "Stochastic Quantization of Real Symmetric Matrix 
Models and Field Theory of Non-Orientable Strings "to appear in the proceedings 
of the Moscow 6-th Quantum Gravity Seminar, June 12-19, 1995 Moscow. 

11 



[10] K. Ito, Proc. Imp. Acad 20(1944)519. 

K. Ito and S. Watanabe, in Stochastic Dif ferential Equations, ed K. Ito (Wiley, 
New York, 1978). 

[11] M. Fukuma, N. Ishibashi, H. Kawai and M. Ninomiya, Nucl. Phys. 427(1994)139. 

[12] J. B. Hartle and S. W. Hawking, Phys. Rev. 28(1983)2960. 

[13] E. Witten, Comm. Math. Phys. 121(1989)351. 

[14] C. Rovelh and L. Smohn, Nucl. Phys. B331(1990)80. 

[15] D. Zwanziger, Nucl. Phys. B192(1981)259. 

L. Baulieu and D. Zwanziger, Nucl. Phys. B193(1981)163. 

[16] N. Nakazawa, Nucl. Phys. B335(1990)546. 

[17] R. Graham, Z. Phys. B26(1977), 397; 
Phys. Lett. A109(1985)209. 

[18] N. Nakazawa, in preparation. 

[19] P. Renteln and L. Smolin, Class. Quantum Grav. 6(1989)275. 

[20] R. Loll, Nucl. Phys. B444( 1995)619. 

[21] N. Nakazawa, Mod. Phys. Lett. 7(1992)841. 

[22] J. Sakamoto, Prog. Theor. Phys. 70(1983)1424. 

[23] H. Rumpf, Phys. Rev. 33(1986)942. 

[24] N. Nakazawa, Mod. Phys. Lett. A5(1990)2407; 
Prog. Theor. Phys. 86(1991)1053. 

[25] J. Greensite Nucl. Phys. B361(1991)729. 



12 



